Constraints on isocurvature models from the WMAP first-year data 
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We investigate the constraints imposed by the first-year WMAP CMB data extended to higher 
multipole by data from ACBAR, BOOMERANG, CBI and the VSA and by the LSS data from 
the 2dF galaxy redshift survey on the possible amplitude of primordial isocurvature modes. A flat 
universe with CDM and A is assumed, and the baryon, CDM (CI), and neutrino density (NID) 
and velocity (NIV) isocurvature modes are considered. Constraints on the allowed isocurvature 
contributions are established from the data for various combinations of the adiabatic mode and one, 
two, and three isocurvature modes, with intermode cross-correlations allowed. Since baryon and 
CDM isocurvature are observationally virtually indistinguishable, these modes are not considered 
separately. We find that when just a single isocurvature mode is added, the present data allows an 
isocurvature fraction as large as 13 ± 6, 7 ± 4, and 13 ± 7 percent for adiabatic plus the CI, NID, 
and NIV modes, respectively. When two isocurvature modes plus the adiabatic mode and cross- 
correlations are allowed, these percentages rise to 47 ± 16, 34 ± 12, and 44 ± 12 for the combinations 
CI-I-NID, CI-I-NIV, and NID-I-NIV, respectively. Finally, when all three isocurvature modes and 
cross-correlations are allowed, the admissible isocurvature fraction rises to 57 ± 9 per cent. The 
sensitivity of the results to the choice of prior probability distribution is examined. 
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I. INTRODUCTION 



Although the current observational data on the cos- 
mic microwave background (CMB) anisotropy and large- 
scale structure (LSS) seem consistent with a spatially 
flat universe with a cosmological constant (A) filled with 
ordinary matter (baryons, leptons, photons, and neutri- 
nos) and a cold dark matter (CDM) component with 
a primordial spectrum of Gaussian adiabatic density 
perturbations described by a monomial power law [i.e., 
P{k) cx fc"'], it is of interest to consider whether mod- 
els having an appreciable isocurvature component can 
also account for the current data. In this paper we con- 
sider the observational constraints on models for which 
in addition to the adiabatic mode either one, two, or 
three isocurvature modes are allowed as well as non- 
vanishing cross-correlations between the allowed modes. 
We consider only the isocurvature modes present if no 
new physics is assumed — that is, the baryon isocurva- 
ture, the CDM isocurvature (CI), the neutrino isocur- 
vature density (NID), and neutrino isocurvature velocity 
(NIV) modes. The baryon and CDM isocurvature modes 
are not considered separately, because they imprint vir- 
tually indistinguishable perturbations on the CMB. Con- 
sequently, without loss of generality, we consider only the 
CI, NID and NIV modes. 

A simple family of models with only adiabatic pertur- 
bations is defined by the cosmological parameters J^a, ^b, 
^cdm (subject to the constraint fl^ + fib + flcdm = 1), 
h, and the reionization optical depth r. These models 
may be compared to the WMAP data [1,2], extended 
to smaller angular scales by including data from the 



ACBAR [3], BOOMERANG [4], CBI [5] and VSA [6] 
experiments and LSS data from the 2dF galaxy redshift 
survey (2dFGRS), [7]. The analysis by the WMAP team 
indicates a good fit to a model of this sort with the 
following best-fit parameter choices: f^A — 0.74, fit, = 
0.043, /i = 0.73, T = 0.15 [8]. 

The object of this paper is to investigate models where 
the requirement of adiabaticity has been relaxed. Adia- 
batic perturbations result when the stress-energy filling 
the universe obeys a single, spatially-uniform equation of 
state. Initially, on surfaces of constant cosmic density, 
the densities of each of the components are also uniform 
and share a common velocity field. This is the simplest, 
but by no means only possible, assumption about the 
character of the primordial perturbations. The simplest, 
single-field models of inflation yield perturbations that 
are exclusively adiabatic, but more complicated pertur- 
bations are possible from multi-field inflationary mod- 
els. Perturbations in the initial equation of state or in 
the relative velocity between the various components are 
called 'isocurvature' because asymptotically, at vanish- 
ing cosmic time, there is no perturbation in the space- 
time curvature, or more precisely this perturbation is of 
higher order. However, as the universe evolves, the vari- 
ous components contributing to the stress-energy evolve 
differently, leading to perturbations in the density, which 
source perturbations in the gravitational potential, whose 
potential wells attract all the components, leading to a 
gravitational instability. One might say that the distinc- 
tion between so-called adiabatic and isocurvature pertur- 
bations is useful only asymptotically at very early times, 
and that almost any perturbation eventually leads to cur- 
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vature perturbations. 

Since the CMB is imprinted at a rather late time in our 
early cosmic history, at z 'i^ 10^, for comparing to the 
data, it suffices to consider a universe containing only 
baryons, CDM, photons, neutrinos and a cosmological 
constant A. This seems to be the minimal physics at the 
relevant energy scales required to account for the present 
observations. There are many possibilities for how the 
perturbations in these components may have been im- 
printed at early times by physics at higher energy scales 
having other relevant degrees of freedom, but for our pur- 
pose, as long as there are no relic particles (other than the 
CDM) left over from such an earlier epoch, the predic- 
tions for the CMB and LSS may be derived entirely from 
the perturbations in the baryon, CDM, photon, and neu- 
trino components. If modes that diverge at early times 
are excluded, there are only five allowed modes: an adi- 
abatic mode (AD), a baryon isocurvature mode (BI), a 
CDM isocurvature mode (CI) [9,10], a neutrino density 
isocurvature mode (NID) and a neutrino velocity isocur- 
vature mode (NIV) [11,12]. (For a review and further 
references see ref. [12]). 

As long as only quadratic observables are considered, 
the perturbations in the five modes just enumerated may 
be completely characterized by a matrix valued power 
spectrum [13,12] 

(^Xiik) X}{k')^=Aum S'{k-k') (1) 

where the indices /, J = 1, 2, 3, 4, 5 label the modes AD, 
BI, CI, NID, NIV, respectively, and the random vari- 
able Xj{k) labels the amplitude of the Ith mode with 
wavenumber k. Because balancing perturbations in the 
CDM against the baryon density leads to a mode that 
hardly evolves at all, we do not consider the baryon 
and CDM isocurvature separately and opt for considering 
only the CDM isocurvature mode. The diagonal elements 
represent auto-correlations whereas the off-diagonal ele- 
ments represent cross-correlations. These must on physi- 
cal grounds be constrained so that Ajj is always positive 
definite. 

A variety of models have been proposed that gener- 
ate mixtures of adiabatic and isocurvature perturbations 
with non-trivial correlations. In an inflationary context, 
any field of mass small compared to the Hubble scale 
during inflation is disordered in a calculable manner, 
with significant power on very large scales. This mech- 
anism is exploited in the curvaton scenarios, where such 
a field is used to set up a variety of initial conditions 
[14,15]. Another interesting possibility arises if inflation 
is driven by two scalar fields [16,17,13]. Recently, it was 
shown that perturbations in the (dominant) scalar field 
which lie along the direction of its motion correspond to 
adiabatic perturbations, while orthogonal perturbations 
correspond to isocurvature perturbations [18]. An ar- 
bitrary trajectory, in general, generates correlated adia- 
batic and isocurvature perturbations with the magnitude 



of the auto- and cross-correlation amplitudes depending 
on the particular model. Other possibilities for generat- 
ing isocurvature perturbations are discussed in [12]. 

It is interesting to ask whether such models may 
be distinguished observationally. We adopt a model- 
indopcndcnt approach, attempting to constrain the am- 
plitude of correlated adiabatic and isocurvature pertur- 
bations directly from the data. A previous investigation 
[19] estimated the quality of the constraints on isocurva- 
ture modes and their correlations that would be obtained 
from the WMAP and PLANCK data based on the pub- 
lished projected sensitivities of these instruments. It was 
shown that when the adiabatic assumption was relaxed, 
the uncertainties on the determination of the cosmolog- 
ical parameters degraded substantially and large admix- 
tures of isocurvature were allowed. Additional polariza- 
tion data, however, was shown to break the degeneracies 
introduced by isocurvature modes, thereby constraining 
non-adiabatic amplitudes to less than 10% and recov- 
ering accurate parameter estimates. In this analysis a 
fiducial model was assumed and estimates on the accu- 
racy of parameter estimation were obtained by examin- 
ing the second partial derivatives of the log likelihood 
about the most likely model. Here we fully explore the 
entire parameter space using Markov Chain Monte Carlo 
(MCMC) methods with the WMAP data, smaller scale 
CMB datasets and LSS data from 2dFGRS. 

Observational constraints on isocurvature perturba- 
tions from cosmological data have been studied previ- 
ously. Earlier work considered the viability of uncor- 
related admixtures of adiabatic and isocurvature modes 
[20] . Prior to WMAP, limits on correlated mixtures of the 
full set of adiabatic and isocurvature modes were investi- 
gated, but over a smaller set of cosmological parameters 
[21]. Recently, constraints on models in which the adia- 
batic mode is correlated with a single isocurvature mode, 
have been obtained using WMAP data [22,23]. 

In this paper we undertake a full analysis within a 
common framework of all combinations of the four al- 
lowed modes. We investigate the constraints on corre- 
lated isocurvature modes that arise from these datasets 
and study the effect that isocurvature modes have on 
parameter extraction. We consider models in which the 
adiabatic and isocurvature modes share a common spec- 
tral index. Although there exist more general classes 
of models that predict independent spectral indices for 
adiabatic and isocurvature modes, the reason for this re- 
striction is that current data, as we shall see, only weakly 
constrains the most general primordial perturbation with 
a common spectral index for pure modes. Allowing inde- 
pendent spectral indices for each pure mode would fur- 
ther degrade parameter extraction due to complex flat 
directions that would arise. Independent spectral indices 
were considered in [23] but only for a single correlated 
isocurvature mode. We have reserved this extension for 
a future study, when such an analysis will become more 
tractable as the quality of data improves. 

In considering the most general primordial perturba- 
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tion we can address one of our main concerns: whether 
it is possible to set stringent bounds on the total isocur- 
vature contribution to current cosmological data. When 
considering specific models at the origin of hybrid initial 
conditions, other groups have tended to restrict them- 
selves to at most one isocurvature mode and its correla- 
tions with the adiabatic mode [14-17]. In general they 
have found that it is unlikely that a considerable frac- 
tion of isocurvature perturbations is allowed. In this pa- 
per we present a comprehensive study of combinations 
of adiabatic and isocurvature modes. We find that we 
essentially reproduce previous constraints if we only con- 
sider a reduced parameter space with one adiabatic and 
one isocurvature mode. If we incrementally increase the 
number of modes, we find that the constraint on the frac- 
tional contribution of isocurvature modes relaxes. Hence 
in this paper we are able to place bounds on a wide range 
of theoretical models, complementing and extending the 
work presented in [24]. 

The paper is organized as follows. Section II details the 
parameterization and prior probability distribution used 
for our analysis. The results are presented in section 
III for the various possible combinations of primordial 
perturbation modes. In section IV the sensitivity of the 
results to various aspects of the prior probability distri- 
bution is investigated. Section V presents some general 
conclusions. 

The analysis in this paper required certain technical 
advances in order to characterize the posterior probabil- 
ity distribution in a computationally efficient and feasi- 
ble manner. Because of the increase in the number of 
undetermined parameters compared to the case of pure 
adiabatic models and the presence of poorly determined 
degenerate directions, it was not possible to carry out this 
analysis with the same technology used for MCMC chains 
for the pure adiabatic models. The techniques used to 
tune the MCMC chains and diagnose their convergence 
have already been described in a separate publication 
[25]. In an Appendix we detail the modifications made 
to the publicly available DASh software [26,27] used for 
calculating theoretical model spectra when isocurvature 
modes and their correlations arc admitted. With this 
software, cosmological models are first evaluated on a 
grid in parameter space. Subsequently, cosmological 
models are evaluated by interpolation over this grid. 

II. PARAMETER SPACE AND DATASETS 

A. Cosmological parameters 

We assume four cosmological fluids: electromagnetic 
radiation p^, baryons pb, non-relativistic cold dark mat- 
ter pc, and massless neutrinos p^, all contributing to the 
total energy density. For the unperturbed model, we de- 
fine oji = Clih^ and fli = STrGpi/SiJo where i labels 
the components, pi indicates the energy densities, and 



Hq = 100 /ikms""'^ Mpc""'^ is the present Hubble con- 
stant. A contribution from a cosmological constant f^A 
is also included. We restrict our analysis to flat uni- 
verses with h > 0.5. We consider the following param- 
eter ranges: 0.03 < tOc < 0.3, 0.01 < Wfc < 0.057, and 

< Oa < 0.9. There is some evidence from a compar- 
ison of adiabatic models with the WMAP data that an 
unconventional ionization history is preferred [8] . In par- 
ticular a large optical depth t w 0.17 seems to indicate 
an epoch of reionization which may be much earlier than 
previously supposed. We therefore allow a wide range of 
admissible optical depth, < r < 0.3. 

B. Mode parameterization 

We characterize the adiabatic and isocurvature pertur- 
bations using the covariance matrix 

(X,(k) X}(k')) = Aij{k)6^{k - k'), (2) 

which provides a complete description of the primordial 
perturbations for Gaussian initial conditions [12]. Here 

1 may take the values AD, CI, NID, NIV, or some sub- 
set thereof containing N pure modes. As previously dis- 
cussed, the baryon isocurvature mode need not be con- 
sidered separately, so we do not include it. We assume 
that the underlying power spectra are smooth over a wide 
range of k, so that the diagonal and off-diagonal elements 
of the matrix A can be parameterized as 

Au{k) = A,jk-'^, = (3) 

in terms of a set of amplitude parameters Aj j and power- 
law shape parameters nj. The cross-correlation spec- 
tral indices are set to the mean of the corresponding 
auto-correlation spectral indices to ensure that A re- 
mains positive semi-definite over a wide range of scales. 
Our convention for the shape parameter is that scale- 
invariant adiabatic and isocurvature spectra correspond 
to n/ = 1. In practice this is achieved by choosing the 
random variable, Xj, that defines the initial power spec- 
trum as Xad = S-y,Xci = k^6c, Xnjd = k^5i,, and 
Xniv = k'^Vv for the different modes, where 5-y,5c and 
5i, are the primordial photon, CDM and neutrino den- 
sity contrasts, respectively, and Vi, is the neutrino veloc- 
ity perturbation, at early times. We assmne a constant 
spectral index ns for all the modes, so that 

Aij{k) = Aij k""' . (4) 

We parameterize so that 

Aij oc zij (5) 

where 

N 

tr{zz^) = J2 = 1 (6) 
/,j=i 
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and the proportionality constant that specifics the over- 
all amplitude is defined below. The matrix elements z/j 
indicate the fractional contributions of the various auto- 
and cross-correlations. There are D = N{N + l)/2 in- 
dependent correlations and the coefficients zjj cover a 
unit sphere of dimension d — D — 1 on which wo use the 
measure corresponding to the usual volume element. 

For the MCMC random walk the coordinates zjj are 
not useful. Wc therefore map the surface of the d- 
dimensional sphere 5"^ onto a d-dimensional ball B''^ using 
the volume preserving mapping 



riO) 



d6 sin 



(d-l) , 



1/d 



(7) 



where 9 is the angular coordinate with respect to the 
north pole of S"', the point corresponding to a pure adia- 
batic model and r is the radial coordinate of B'^. Conse- 
quently, a random walk in B"^ (subject to the constraint 
6 < tt) may be lifted to S'^ by the inverse mapping. More 
explicitly, 

sin 6* Wi sm9 ww-i 



-222 



: Smt' , . . . , Z2N 

r 



znn 



SmO W2N-2 



V2 



smf 



5 

r 



(8) 



and zij = zji, where Wi are the Cartesian coordinates of 
the Euclidean space into which B"^ had been embedded, 
r = |w|, and 6 = 9{w). In practice, we choose / = 1 to 
label the adiabatic mode, with 2. . .N labeling the isocur- 
vature modes. This parameterization covers the space of 
all symmetric N x N matrices. We implement positive 
definiteness by assigning a zero prior probability to any 
proposed matrix z having a negative eigenvalue. Since 
the eigenvalues can be computed quickly, any proposed 
steps in the MCMC sequence are efficiently rejected be- 
fore the corresponding cosmological model is calculated. 

The modes are first normalized to give equal power to 
the CMB anisotropy summed from £ = 2 through £ = 
2000 inclusive. The power of an auto-correlated mode / 
is defined as 



2000 



J2{2£+1)C'/. 



(9) 



t.=2 



For the cross-correlations the geometric means of the re- 
spective renormalization factors is used. This defines the 
fractional contributions zjj in terms of the physically ob- 
servable power in each mode. 

We now define the constant of proportionality in oqn. 5. 
We choose to sample the total power in the CMB 
anisotropy, A = X)^(2£-|-l) Q, where Cg = ^ijC^'^ , 
since it is one of the best determined quantities from the 
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FIG. 1. Shape of prior distribution for the observ- 

ables 2(1.1) and Z{ij) ■ Our prior distribution of matri- 
ces subject to the constraint tliat tr(zz^) — 1 and that 
Z{l^j) have only positive eigenvalues, in the absence of data, 
would give these posterior distributions for the fractional 
auto-correlations i) (dashed) and cross-correlations 2:(i.,i> 
(dot-dashed) for A'' = 4. The solid curve shows the distri- 
bution for Z(i^Yi that would result if the positive definiteness 
requirement were suppressed. 



CMB data and hardly varies among models lying in the 
high-likelihood region of parameter space. To do so we 
perform a second renormalization so that 



Aij = A 



^^(2£ + l) zjjCl' 



ZlJ 



(10) 



Except for the second renormalization, the prior dis- 
tribution of Z is invariant under the transformations 

Z^OZO-^ (11) 

where O £ SO{N). Other parameterizations may be con- 
templated. For example, Z may be expressed in the form 

(12) 



Z = ODO 



-1 



where O is orthogonal and D is diagonal. This is es- 
sentially the approach proposed in [21].^ The space of 
orthogonal matrices has a natural volume element, be- 
cause the requirement that volume be invariant under 
left (or right) translation fixes the form of the volume el- 
ement up to a constant factor. It is, however, less obvious 
what measure should be chosen for the space of diagonal 
matrices. 

In our approach the effective measure for D is propor- 
tional to 



^Actually, in [21] a Cartesian measure is chosen for where 
O = exp[QiEi]. This measure, however, is not invariant under 
group translations. 
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1 ni^^-^-^ 



(13) 



KJ 



[28] where A/ labels an eigenvalue of D. The second 
factor in the above expression gives rise to the well know 
phenomenon of "level repulsion." Eigenvalues of a ran- 
dom matrix do not follow Poisson statistics but rather 
are anti-bunched. This level repulsion greatly disfavors 
models for which one eigenvalue dominates, because in 
this case all the others arc pushed toward zero. In sec- 
tion IV we shall consider rescalings to counteract this 
phenomenon. The effect of this level repulsion is illus- 
trated in Fig. 1 where the prior distributions for zij are 
plotted for N = 4. The solid line corresponds to p{zij) 
for all modes if no positive definite requirement is placed 
on the matrix. For smaller N the effect is less extreme. A 
further phase-space effect arises from the positive semi- 
definite requirement, which constrains Z to have positive 
eigenvalues. This results in a phase space suppression 
at low auto-correlation amplitudes zn and high cross- 
correlation amplitudes zjj,I ^ J. The resulting prior 
distributions for these modes are shown in Fig. 1. 

Using these conventions we define the non-adiabatic 
fraction as 



/iso = 



ziso 



ZlSO 



(14) 



^(AD,AD) 

with the isocurvature contribution to the data given by 



^ISO 



'(AD,AD) ' 



(15) 



which includes the contribution from cross-correlation 
modes. 

Another possibility, in terms of amplitudes, is to con- 
sider v^/iso as an alternative candidate for the isocurva- 
ture fraction. We note that other measures of the non- 
adiabatic fraction are possible as well, particularly in the 
presence of several modes and their correlations. 

Because of the possibility of cross-correlations, the to- 
tal power is not simply the sum of the adiabatic and the 
isocurvature powers. It is the amplitudes that add, the 
power being the square of their sums. Consequently, con- 
structive or destructive interference or a combination of 
these may take place. It is therefore theoretically possi- 
ble to have virtually identical adiabatic and total power 
in the presence of a large isocurvature contribution. 



C. Datasets 

For the CMB data, we use the WMAP first-year Cf 
and Cf^ data [1] covering 2 < £ < 800 using the like- 
lihood function in [31] and extend the data for Cf to 
£ > 800 from the ground-based and balloon-borne ex- 
periments ACBAR [3], BOOMERANG [4], CBI [5], and 
VSA [6], using the compilation in [29,30] with the covari- 
ance matrix for £ < 800 excised. 
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FIG. 2. Pure adiabatic model. Posterior distributions 
for the cosmological parameters using CMB data (dot-dashed) 
and CMB-I-LSS data (dashed). 



For the large-scale structure data, we use the galaxy 
power spectrum Vg(k) as measured in redshift space by 
the 2dFGRS experiment [7] at an effective survey redshift 
of 2; = 0.17. We use only data in the linear regime 0.01 < 
k/{hMY>c-^) < 0.15 and assume that the distortion due 
to redshift distortion and the bias are independent of 
scale. Following [31] we rescale according to 



Vg{k)U=o.i7 = f{b,(3) Pm{k)U=o 



where the rescaling factor 



/(6,/3) = 6^(l + ^/3e// + J/3e'//) 



(16) 



(17) 



is given in terms of the redshift distortion factor /3 and 
bias b. We approximate the bias as 6 » Vl^-;^ / (3 and take 
P^f^ = 0.85/3 as the correction to /? at redshift z = 0.17. 
For the LSS data from 2dFGRS data (3 is an additional 
parameter with a Gaussian prior 0.43 ±0.077 [32] broad- 
ened by 10% as in [31] to allow for possible error in /3e//- 



III. RESULTS 

In this section wc present our results, first consider- 
ing the pure adiabatic model, essentially repeating the 
analysis of the WMAP team, and then systematically 
introducing isocurvature modes and their correlations. 
We first consider the addition of a single isocurvature 
mode, then the addition of pairs of isocurvature modes, 
and finally the addition of all three allowed isocurvature 
modes. 

Fig. 2 shows the marginalized posterior probability dis- 
tributions for the cosmological parameters, both for the 
CMB dataset and for the combined CMB-I-LSS dataset, 
for pure adiabatic models. The median values and 68% 
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CMB 



CMB+LSS 



AD+CI 



AD+NID 



AD+NIV 



UJb 
U)c 

Oa 

Us 

T 

P 

h 
b 



0.024 ± 0.001 0.023 ±0.001 
0.13 ±0.01 0.120 ±0.006 
0.71 ±0.03 
0.97 ±0.03 

'-'•-'■^-0.06 

0.50 ±0.05 



^•0^-0.08 

0.99 ±0.03 
0.15 ±0.07 



n o-i_|_0.08 

0.69 ±0.05 



0.29 ±0.03 
0.70 ± 0.03 
0.95 ±0.08 



TABLE I. Pure adiabatic model. Median parameter 

values and 68% confidence intervals using CMB data (column 
1) and CMB+LSS data (column 2). 
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FIG. 3. Adiabatic plus CI. We indicate the marginalized 
posterior distributions of the various parameters for the adi- 
abatic (AD) and CDM isocurvature (CI) mixed models with 

correlations allowed (solid) using CMB+LSS, with the result 
for the pure adiabatic model (dashed) included for compari- 
son. (i,j) denotes j> . 



confidence intervals for each parameter are given in Ta- 
ble I. Our results are consistent with the WMAP anal- 
ysis [8]. The density parameter values Ui inferred from 
CMB and LSS measurements are consistent with those 
inferred using other techniques, such as measurements 
of the baryon density from observations of primordial 
light element abundances [33] and of the energy density 
of a nc\a,ativc pressure component inferred from Type la 
supernovae [34]. Although the data suggests an earlier 
reionization indicated by a large optical depth r, this pa- 
rameter is not very well constrained by the CMB. The 
tilt of the adiabatic power spectrum is also well mea- 
sured, with a nearly scale-invariant spectrum favoured. 
Table I also lists several derived parameters such as the 
matter fraction O^, the dimensionless Hubble parameter 
h, and the bias b. 



Da 

Us 
T 

Z(Cl,Cl) 

^(NID.NID) 

•2(NIV,NIV> 

2(AD,CI) 

^(AD,NID> 

2^(AD,NIV> 

/iso 

h 
b 



0.026 ±0.003 
0.116 ±0.007 
0.74 ±0.04 
1.04 ±0.06 

0.46 ±0.05 



qq+0.009 
0.10 ±0.06 



-0.08 ±0.05 



0.15 ±0.08 
0.13 ±0.06 

0.26 ±0.04 
0.74 ±0.05 
0.96 ±0.08 



0.025 ±0.002 
0.118 ±0.007 
0.73 ±0.04 
1.02 ±0.05 

n u+o-08 

U.i4_o.06 

0.48 ± 0.06 



n 007+0.003 
u.yy (_o,oo6 



0.05 ± 0.04 



-0.03 ±0.03 



0.08 ± 0.05 
0.07 ±0.04 

0.27 ±0.04 

0.73 ±0.05 
0.94 ±0.08 



0.026 ±0.002 
0.112 ±0.007 
0.71 ±0.03 
1.02 ± 0.04 
21+0 06 

0.47 ±0.05 



0.99 



+0.008 
-0.02 



0.07 



+0.06 
-0.04 



0.08 ±0.08 

0.15±8:J? 
0.13 ±0.07 

0.29 ±0.03 

0.69 ±0.04 
1.00 ±0.10 



TABLE II. Adiabatic plus one isocurvature mode. 
Median parameter values and 68% confidence intervals for 
mixed models with the adiabatic mode (AD) and a sin- 
gle isocurvature mode (columns 1-3), with the CMB+LSS 
dataset used throughout. 



A. Adiabatic mode correlated with single 
isocurvature modes 

We present results for each of the cold dark matter 
isocurvature (CI), neutrino isocurvature density (NID) 
and neutrino isocurvature velocity (NIV) modes corre- 
lated with the adiabatic mode. The median values and 
the 68% confidence intervals of the cosmological param- 
eters and mode amplitudes are given in Table 11 for each 
of the above mode combinations. Of immediate inter- 
est arc the relative powers that are permitted for the 
pure isocurvature modes and their cross-correlations. In 
all three cases the preferred models are dominated by 
the adiabatic mode, with a small non-adiabatic fraction, 
/iso ^ 10%, tolerated. 

The relative powers z/j of the individual isocurvature 
auto-correlation and cross-correlation modes are consis- 
tent with zero at the Icr to 2(j level. The marginalized 
parameter distributions for the CI, NID and NIV modes 
are shown in Figs. 3, 4, and 5, respectively. The best fit 
mixed models (which have a 10% non-adiabatic contri- 
bution) do not fit the data significantly better than the 
best-fit adiabatic models. An illustration of the flat di- 
rection that arises in correlated adiabatic and CI models 
was given in [18]. 

Including a single correlated isocurvature mode does 
not significantly change the cosmological parameters 
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FIG. 4. Adiabatic plus NID. Distributions as described 
in Fig. 3, for correlated adiabatic and neutrino density isocur- 
vature models. 
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FIG. 5. Adiabatic plus NIV. Distributions as described 
in Fig. 3, for correlated adiabatic and neutrino velocity isocur- 
vature models. 



from those in pure adiabatic models, as shown in Table 
II. The rig distribution is slightly broadened and shifted 

to larger values, and ojf, is slightly increased. For the NIV 
mode T shifts toward larger values, and the distribution 
appears cut-off by the limit at t = 0.3 from the prior. 
This is because the NIV mode adds most of its power on 
small scales 1°) and its contribution to the CMB is 
damped by the larger r. 



B. Adiabatic mode correlated with two isocurvature 
modes 

We now study the effect of introducing two isocurva- 
ture modes correlated with the adiabatic mode and each 
other. The three possible combinations are CI+NID, 
CI+NIV, and NID+NIV. In all cases the isocurvature 



^iso 
/iso 

h 
b 



0.66±Si« 



o.46±S:}« 



0.47 ±0.16 0.34 ±0.12 0.44 ±0.12 



0.26 ±0.03 
0.78 ±0.06 
1.00 ±0.10 



0.27 ±0.03 
0.72 ±0.05 
1.13 ±0.14 



0.27 ±0.03 
0.75 ± 0.04 
1.21 ±0.16 



TABLE III. Adiabatic plus two isocurvature modes. 

The same conventions as in Tabic II arc used, for mixed mod- 
els including two correlated isocurvature modes. 
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FIG. 6. Adiabatic plus CI+NID. Distributions as de- 
scribed in Fig. 3 for the adiabatic mode correlated with both 
the CDM and neutrino density isocurvature modes. 




FIG. 7. Adiabatic plus CI+NIV. Distributions as de- 
scribed in Fig. 3 for the adiabatic mode correlated with both 
the CDM density and neutrino velocity isocurvature modes. 



FIG. 8. Adiabatic plus NID+NIV. Distributions as de- 
scribed in Fig. 3 for the adiabatic mode correlated with both 
the neutrino density and velocity isocurvature modes. 



fraction increases significantly, to /iso ~ 47%, 34% and 
44% respectively, as can be seen from Table III. Tlie 
posterior distributions are shown in Figs. 6, 7 and 8. 

The isocurvature contribution is greatest for the 
CI+NID models, although apart from an increase in 
the baryon density, the cosmological parameters are not 
shifted significantly from their median values in the adi- 
abatic case. Despite the large isocurvature fraction, the 
best-fit models have both positive and negative mode 
correlations which cancel to give a small non-adiabatic 
power, as illustrated in Fig. 9, leaving the adiabatic mode 
with a power closely matched to the total spectrum. To 
gain insight into these non-adiabatic cancellations, we 
utilise the covariance matrix of the distribution to per- 
form a principal component analysis about a high likeli- 
hood mixed model with parameters Xi. We calculate the 
eigenvalues, Aq,, and eigenvectors, y^, of the covariance 
matrix, 

dj = {hxj) - {xi){xj) (18) 

where Xi = xi/ Xi are normalized by their fiducial values. 
The flattest direction, yo, along which changes in the 
CMB and matter power spectra are negligible, has per- 
centage error \Syo\ = 100 %/v^ = 90%. Fig. 10 shows 
the cancellation in the CMB temperature spectrum in 
this direction, dominated by a subset of seven parameters 
that comprises the spectral index, rig, and the power in 
the six mode correlations. The increased auto-correlation 
power on large-scales that results from a low spectral in- 
dex is compensated almost exactly by the negative cross- 



correlation spectra (ad,ci) and (ad.nid). This highlights 
the important role played by these cross-correlations in 
allowing a large isocurvature contribution. 

The Cl-f NIV models allow a lower isocurvature frac- 
tion than those with CI-I-NID, but the baryon density, 
spectral index and optical depth r are further increased 
and the redshift distortion parameter reduced from the 
pure adiabatic case. As observed in the previous section, 
the increase in r is correlated with an increase in ampli- 
tude of the NIV mode. The increased NIV contribution 
is also strongly correlated with an increase in the values 
of the spectral index, n^, and the baryon fraction, w^, 
which are all associated with a flat direction in parame- 
ter space, as we will show in the next section. 

The cosmological parameters change more signiflcantly 
from their adiabatic values when both neutrino isocur- 
vature modes are introduced, shown in Fig. 8. The 
baryon density, spectral index and reionization param- 
eter are significantly increased, while the redshift distor- 
tion parameter is further reduced. The best-fit models 
have positive- valued cross-correlations, adding construc- 
tively as indicated in Fig. 9 to produce a significant non- 
adiabatic amplitude. The lower value of /? in these mod- 
els is then consistent with the reduced adiabatic power 
and higher normalization factor f{b, (3) required to fit the 
LSS data. 



8 



6000 F 




FIG. 9. CMB temperature spectra for best fit mixed 
models. (Top) The CMB angular power spectrum, Cj'^ , for 
the best-fit mixed model (solid) including AD+CI+NID 
modes with parameters (wj,, Wd, Qa, ris, r, /iso) = (0.033, 
0.13, 0.75, 0.99, 0.16, 0.62). The adiabatic (dot-dashed) and 
non-adiabatic (dashed) contributions arc included. (Bottom) 
The same for models with AD-|-NID-|-NIV, for parameters 
(0.037, 0.11, 0.74, 1.11, 0.22, 0.56). The CMB data is in- 
cluded. 



C. Adiabatic mode correlated with all three 
isocurvature modes 

We now consider a general perturbation comprising a 
linear combination of the AD, CI, NID and NIV modes. 
Table IV gives statistics for the cosmological parameters 
and relative mode amplitudes for both the CMB and 
combined CMB-I-LSS datasets, and these distributions 
are plotted in Fig. 11. 

Wc observe that the adiabatic mode is no longer com- 
pletely dominant. The general model includes a mixture 
of adiabatic and isocurvature modes with /iso = 60% for 
CMB data alone which decreases slightly to /iso = 57% 
for the CMB-I-LSS dataset. The mixed model family 
include models with baryon fractions greater than that 
inferred from nucleosynthesis measurements as well as 
large optical depths and large departures from scale- 
invariancc. These departures from the pure adiabatic pa- 
rameter values are all correlated with larger isocurvature 
fractions and have broad distributions. These shifts are 
associated with a degenerate direction in parameter space 
which we investigate below. The remaining cosmological 
parameters Wc and are consistent with their values 
in the adiabatic case for the CMB-I-LSS dataset. This is 
also true for the derived parameters, though there is a 
preference for slightly smaller values of Q.m and slightly 
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FIG. 10. Derivative along flattest direction of the 
CI-|-NID distribution (dot-dashed) is decomposed into 
variations in the mode contributions, (l,J), and the overall 
spectral index ris- The shaded region corresponds to the un- 
certainty of the WMAP measurement including cosmic vari- 
ance. 
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0.19 + 0.11 
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0.26 + 0.03 
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TABLE IV. Adiabatic plus all three isocurvature 
modes. The same conventions are used as in Table II, 
for mixed models with all three correlated isocurvature 
modes (AD+ISO). Results are included for both CMB and 
CMB+LSS datasets, and for pure adiabatic models (AD) for 
comparison. 

larger values of h. 

In Fig. 12 we plot the CMB and matter power spec- 
tra for an adiabatic model and a mixed model with high 
likelihood, with parameters (w^, cjc, Oa, rig, r, /3, /iso) 
=(0.041, 0.13, 0.75, 1.06, 0.28, 0.37, 0.60). The adiabatic 
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and mixed model spectra are essentially indistinguish- 
able, despite the adiabatic mode contributing only 40% 
of the total rms power of the mixed model spectra. 

Although the posterior distributions for mixed mod- 
els and pure adiabatic models hardly overlap, the high- 
likelihood models are connected in parameter space. 
That is to say, the goodness-of-fit, x^, changes almost 
monotonically over a variation of parameters that in- 
terpolates linearly between these models, as shown in 
Fig. 13. This means that the posterior distribution is 
not bimodal but rather flat between the high likelihood 
adiabatic model and the high likelihood mixed model. 
A principal component analysis of the covariance matrix 
of the distribution, calculated about the high likelihood 
mixed model shown in Fig. 12, indicates that the flattest 
direction yo has an uncertainty greater than 100%, with 
three eigen-directions having uncertainties greater than 
50%. 

To elucidate what physical effects are responsible for 
these flat directions, we search for degenerate directions 
that arise from smaller subsets of the full parameter set. 
We have identified two such flat directions. The flrst is 
the degenerate direction described in Section III A in- 
volving seven parameters, which is again present here. 
This flat direction accounts for the negative correla- 
tions, 2:(AD,ci> and 2;(ci,NiD) , and the positive correlation 
•2(AD,NiD) • By performing a principal component analy- 
sis of these seven parameters, we find a very similar flat 
direction to the one shown in Fig. 10. 

The second flat direction, that involves four param- 
eters, is associated with a large shift in the value of 
the baryon density. By studying the two-dimensional 
marginalised distributions shown in Fig. 14, it is clear 
that the spectral index and the adiabatic and neutrino 
isocurvature velocity mode contributions are strongly 
correlated with the baryon density. Performing a princi- 
pal component analysis of the covariance matrix of this 
four-parameter set, reveals a flat direction with an uncer- 
tainty of 50%. This cancellation is illustrated in Fig. 15, 
where the baryon density is increased by 8%, leading to 
the suppression (boosting) of even (odd)-numbered adi- 
abatic acoustic peaks. This shift is compensated by in- 
creasing the power in the NIV mode by 29%, reducing 
that in the AD mode by 13% and increasing the overall 
spectral index by 2%, leaving the CMB temperature 
spectra unchanged to within the WMAP experimental 
error bars. It is not surprising that the neutrino velocity 
mode provides the oscillatory cancellation as the baryon 
density increases, given that it is approximately in phase 
with the adiabatic mode. 

We now consider the statistical significance of these 
results. The mixed model reduces the total by 5. Un- 
der the approximation of a linear parametric dependence 
(which in this case is probably not a very good approxi- 
mation), one would expect the reduction in to be gov- 
erned by a distribution with nine degrees of freedom 
(i.e., equal to the number of extra parameters). Conse- 
quently, a reduction by only 5 where simply fitting the 



!l /' 


u 
\\ 


M /; 
' ' /' 


\\ 
\\ 








' ' /' \ 


\ 


! // ^ 
! A 


\\ 


ij v 





0.08 0.12 


0.1 6 






' V 




1 /\ 

/ / 






\ 








0.8 1.0 1.2 1.4 0.0 0.1 0.2 0.3 0.2 0.4 0.6 




/ \ 


<c/,c/> - 


- / \ 





0.0 0.7 0.4 0.6 0.8 1.0 0.0 0.0 0.4 0.6 





^<NID,NID> - 




\\ 


1/ 


\\ 




\\ 




\\ 


-J 






CO 0.2 G.4 0.6 0.0 0.2 0.4 0.6 0.8 

<AD,a> ' ' n 1<AD.NID> ' n t<AD,NIV> 





<NID,NIV> 






V, 



-0.4 -0.2 0.0 0.2 0.4 -0,4 -0.2 0.0 0.2 0.4 



-0.4-0.2 0.0 0.2 0.4 



FIG. 11. Posterior distributions for Adiabatic plus 
CI-|-NID-|-NIV. The distributions for models with corre- 
lated adiabatic and three isocurvature modes are shown, using 

the CMB only (dashed) and (CMB+LSS) (sohd). Where rele- 
vant the posterior for the pure adiabatic models (dot-dashed) 
is shown for comparison. 



noise would lead to a reduction by 9, docs not constitute 
evidence in favor of adding the extra parameters. 

We conclude that although our results demonstrate 
that the data allows for large amounts of isocurvature, 
we have not found any evidence indicating that the mod- 
els with more parameters offer a statistically significantly 
better fit. The same analysis may be applied to the ad- 
ditions of a single mode and pairs of modes. Adding a 
single mode adds two parameters, and we find a reduc- 
tion of the total by 1 for the NIV mode, and less than 
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FIG. 13. foi" models connecting adiabatic and 
mixed models, calculated along a path connecting the 
best-fit adiabatic model xi with the best-fit mixed model X2, 
calculated at position x = Ax2 -|- (1 — A)xi. 
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FIG. 14. Two dimensional marginalized distribu- 
tions. Correlations between the baryon density ujb, the spec- 
tral index Ua and the contributions of the pure adiabatic and 
neutrino velocity isocurvature modes are illustrated. de- 
notes J) . 



FIG. 12. Power spectra of best-fit mixed models 
for three isocurvature modes. The angular power spec- 
trum of CMB anisotropics, Cj, the temperature-polarization 
cross correlation spectrum, Cj^ , and the galaxy power spec- 
trum, P{k), for the best-fit adiabatic model (dashed line) and 
best-fit mixed model with /iso = 60% (solid line), with the 
auto-correlation contributions to the mixed model included. 
The CMB data points arc overplotted and the sizes of the 
2dFGRS error bars are indicated. 



1 for the CI and NID modes. In none of these cases is the 
improvement in fit statistically significant. Similarly, for 
pairs of modes, five parameters are added, and the total 
decreases by approximately 1, 1, and 2 for the combi- 
nations CH-NID, CI-FNIV, and NID-hNIV, respectively. 

It is possible in principle that a statistically significant 
improvement in fit might be obtained if the exponents 
in the power laws for the various modes were allowed 
to vary independently and/or additional effects such as 
tensor modes were allowed. We did not investigate this 
possibility. 



IV. SENSITIVITY OF POSTERIOR 
DISTRIBUTIONS TO CONSTRAINTS AND 
PRIOR DISTRIBUTIONS. 

The results of the previous section were obtained using 
Bayes' theorem 

Pposierior{(^par\^data^ L(pCdata\o^par^ Ppriori^^^par) (19) 

where apar spans the parameter space of models and 
L{Tidata\ctpar) IS the likclihood of the model labeled by 
Cipar given the CMB (or CMB+LSS) data indicated by 
Xdataj subject to the prior distribution over the space of 
models Pprior{otpar)- The posterior distribution is com- 
puted using MCMC techniques outlined in [25] and ref- 
erences therein. Certain prior distributions in the form of 
constraints, such as the positivity of various parameters 
that physically cannot be negative or the positive def- 
initeness of the matrix-valued power spectrum, do not 
merit discussion. For most parameters uniform priors 
were assumed, following customary practice. 

In this section we investigate the sensitivity of our 
conclusions (i.e., the posterior probability distribution) 
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FIG. 15. Degenerate direction in ujt, Us, ^{AD.AD); 
and ^(Niv,Niv> The total derivative along the flattest de- 
generate direction (dot-dashed) is decomposed into variations 

in the four parameters above (all sohd), which cancel to within 
the limits allowed by WMAP error bars (shaded). 

for those variables for which the choice of prior was im- 
portant. In particular, we investigate the sensitivity to: 
(1) The allowed range for the reionization optical depth 
r, (2) the priors for the rcdshift distortion parameter 
f3 — and the bias b, (3) the prior for ujh, that is 

whether a uniform prior is employed or whether informa- 
tion from other independent determinations of u)b based 
on nucleosynthesis are employed, and (4) the choice of 
parameterization (or measure) for Aij{k). In all cases 
we consider the = 4 example, including all possible 
isocurvature modes. 



A. Alternative Priors for Cosmological Parameters 

1. Constraint on optical depth 

Following the WMAP team analysis [8] of pure adia- 
batic models, we limited the reionization optical depth to 
the range < r < 0.3. When isocurvature modes were 
included, in particular the NIV mode, we observed that 
the distribution for t seemed to concentrate toward the 
endpoint at r = 0.3. To investigate the sensitivity to the 
value of this cutoff, we considered an alternative upper 
limit at T = 0.5. Fig. 16 compares the posterior distribu- 
tions for these two upper bounds, showing those param- 
eters for which the difference is greatest. Even with the 
upper bound raised to r = 0.5 the posterior distribution 
still accumulates there. However, the effect on other pa- 
rameters is mild, with only a slightly larger isocurvature 
contribution at /iso = 0.60 ± 0.10 permitted. 

2. Priors from galaxy redshift survey data 

We have included data on the galaxy power spec- 
trum from the 2dFGRS survey, incorporating a Gaus- 
sian prior for the redshift distortion parameter with 
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FIG. 16. Effect of changing the upper limit for r 
from 0.3 (solid) to 0.5 (dashed). Only the posterior dis- 
tributions for parameters where the difference is significant 

arc shown; the distributions for the remaining parameters are 
unchanged from their r < 0.3 values. 

[3 = 0.43 ± 0.077, that is derived from a measurement 
of the redshift distortion effect observed in the galaxy 
survey [32,31]. The WMAP team [8] have argued for an 
independent prior on the bias of 6 = 1.04±0.11, obtained 
by computing the bispectrum of the 2dFGRS survey [35]. 
In the previous section we did not impose this prior on the 
bias since it was obtained under the assumption of Gaus- 
sian initial fluctuations. Furthermore, an implicit prior 
is placed on (and therefore JIa) when it is combined 
with the prior on fi. We found a higher value for the bias, 
b = 1.3 ± 0.2, obtained from constraining three isocurva- 
ture modes with the CMB-I-LSS dataset. A larger bias 
generally allows for a larger isocurvature fraction, so we 
now include the additional determination of b to deter- 
mine its effect on the allowed isocurvature contributions. 
The effect is small, with parameter value distributions 
virtually the same and the isocurvature contribution de- 
creasing only slightly from /iso ~ 0.57 to 0.53. 

We have also investigated the role of the prior on /3 in 
constraining the ratio of the galaxy power spectrum to 
the matter power spectrum. f{b, /3). By relaxing the prior 
on P, we allow / to vary freely so that only information on 
the shape of the power spectrum is utilised. The resulting 
parameter distributions are shown in Fig. 17, where / has 
been optimised at each step in the chain. We observe that 
larger isocurvature fractions, /iso = 0.65 ± 0.09 are al- 
lowed. In particular the NIV mode contribution increases 
to 2(Niv,Niv> = 0.40 ±0.15. We have shown that such an 
increase is correlated with a reduction in the amplitude 
of the adiabatic mode. Such models were previously dis- 
favoured by the prior on (3 because of their resulting low 
matter power spcctnmi amplitudes, but are permitted 
when the normalisation of the matter power spectrum 
is free to vary. More models along the flat NIV direc- 
tion are explored, thereby permitting larger values of Hs 
and cub. We also observe weaker constraints on iVc and 
This test demonstrates that the redshift distortion 
measurement is very informative. 
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FIG. 17. Effect of discarding LSS amplitude infor- 
mation. Marginalized parameter distributions obtained us- 
ing only the shape information from the LSS data (dashed) 
compared to those distributions obtained using the prior 
(3 = 0.43 ± 0.077 (solid). 
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-0.06 ±0.11 
0.0008 ±0.11 

0.06 ± 0.09 
-0.10 ±0.09 
-0.003 ± 0.04 

0.04 ± 0.05 

0.41 ±0.14 
0.31 ± 0.09 



0.61 ±0.15 
0.23 ±0.11 

0.30 ±0.12 

+ 0.14 
-0.11 
+ 0.12 
-0.10 



0.28 
-0.12 



+ 0.10 
- 0.12 



0.11 : 

0.19 ±0.11 
-0.18 ±0.10 
-0.09 ± 0.08 

0.16 ±0.08 

0.57 ±0.09 



3. Prior on baryon density 

In the previous section, when three possibly corre- 
lated isocurvature modes were allowed, we found Ub = 
0.041 ± 0.006 for the baryon density, a value discrepant 
with nucleosythesis based determinations of ujt, which in- 
dicate Ub = 0.022±0.002 [33]. We investigate the effect of 
replacing the uniform prior on LUb with a prior that incor- 
porates this complementary information. Fig. 18 shows 
the posterior distributions with and without the nucle- 
osythesis determination. The median values and 68% 
confidence levels are shown in Table V. The posterior for 
the baryon density overlaps with the prior nucleosythesis 
distribution but values as large as Ub = 0.030 are still 
permitted. There is a large reduction in the isocurvature 
contribution, in particular the NIV mode and its correla- 
tions, with the overall non-adiabatic fraction reduced to 
a median value of /iso = 0.31. The cosmological param- 
eters correlated with the NIV amplitude, such as n^, r 
and /3 shift back toward their pure adiabatic values. 

B. Dependence on mode parameterization 

We examine the sensitivity of the posterior distribution 
to the choice of prior distribution for the coefficients of 
the matrix Aij = Ajj fc"^. As detailed in section IIB, af- 
ter normalizing the modes according to their contribution 
to the total CMB power from £ = 2 through £ = 2000, 
we parameterized Ajj as Ajj oc zjj where the constraint 
II 2; II = 1 was imposed and the uniform spherical measure 
on the resulting sphere was assumed. A uniform prior 
for the constant of proportionality was assumed. 



ilm 0.29 ±0.03 0.26 ±0.03 
h 0.70 ± 0.04 0.80 ± 0.05 
b 1.1 ±0.1 1.3 ±0.2 

TABLE V. Effect of incorporating the nucleosynthe- 
sis determination of ojb- Median parameter values and 68% 
confidence intervals as in Table II for mixed models with a nu- 
cleosynthesis prior (BBN, 1st column), compared to the stan- 
dard results with a broad uniform prior (2nd column). The 
CMB±LSS dataset is used throughout. 




0.02 0.04 0.06 0.08 0.12 0.16 0.20 0.60 0.70 0.80 0.8 1.0 1.2 1.4 




-0.2 0.0 0.2 0.4 -0.5-0.3-0.1 0.1 0.3 -0.5 -0.3 -0.1 0.1 0.^0.4 -0.2 0,0 0.2 0.4 



FIG. 18. Effect of incorporating the nucleosynthesis 

determination of iut. The curves show the posterior dis- 
tributions from the CMB±LSS data with a flat prior for uib 
(shown in solid) and a Gaussian prior reflecting the determi- 
nation UJb = 0.022 ± 0.002 (shown as dashed). 
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FIG. 19. Alternative prior uniform in the coeffi- 
cients Aij. The posterior distributions with this new prior 
(dashed) are compared to those of the 'standard' prior de- 
scribed in Section II (solid). 

We now consider directly sampling a uniform distribu- 
tion in the coefficients Ajj. This alternative prior makes 
a large difference for the case of three isocurvature modes, 
for which the allowed isocurvature rises from /iso = 0.57 
with the old prior to /iso = 0.81 with the new prior. 
The posterior distributions are modified significantly, as 
indicated in Fig. 19. 

We can understand this effect by observing that the 
constant of proportionality in eqn. 5, which we now la- 
bel A', can be written as AjjY^^ ■ Under this new 
parameterization where we sample uniformly in Ajj, we 
therefore replace the uniform prior on A' with a prior so 
that P(A') cx (A')" where a = N{N + l)/2 - 1. This 
favours models with high values for A' , corresponding to 
those with large positive and negative isocurvature con- 
tributions Ajj, which may cancel out to fit the data. As 
discussed in the previous section, such models may be 
formed by moving along a degenerate direction consisting 
of the auto- and cross-correlations formed by the CI and 
NID modes with the adiabatic mode, and the spectral 
index Us- A similar direction was illustrated in Fig. 10. 
These models account for the increased contribution of 
both the CI and NID mode auto- and cross-correlations, 
and hence the reduction of the relative AD and NIV mode 
contributions. While the likelihood of such models with 
very high isocurvature is relatively low, the phase-space 
effect due to the modified prior boosts their posterior 
probability. 

Both of the previous priors disfavor models where any 
one mode dominates over the others. This is because of 

the product factor in the distribution given in eqn. (13). 
For the prior presented in section IIB, taking the limit 
-2(AD,AD) 1 causes one of the eigenvalues of zjj to ap- 
proach one while the other (A^ — 1) eigenvalues approach 



FIG. 20. Effect of reweighting to correct for bias 
against pure models. The left panel shows what poste- 
rior for the adiabatic mode contribution 2{AD,AD) would be 
in the absence of data. In the center panel we show how the 
posterior (from the prior in section IIB) (dashed) shifts (to 
the solid) curve) from a re-weighting that would give a flat 
posterior in the absence of data. The right panel shows the 
corresponding shift in uJb resulting from the same reweighting. 



zero. More quantitatively, if ^J(ad,ad) — 

l-e2 the (N-l) 

small eigenvalues cluster within e of zero, giving a factor 
proportional to e^^^~^^^^ in the product in eqn. (13). 
This gives a prior density for -2(ad,ad)5 marginalized 
with respect to the other parameters, proportional to 
gAf(Af-i)/2-i j-j-^ ^Yic neighborhood of Z(ad,ad)- Conse- 
quently, even if the data indicated that all models where 
equally likely — that is, if the data was non-informative — 
our posterior would have a zero of the same order at the 
endpoint Z(ad.ad) = 1' 

We may correct for such bias, or estimate its effect, by 
rescaling the prior so that in the absence of data the pos- 
terior in 2;(AD,AD) would be flat, as shown in Fig. 20. The 
left panel shows the posterior for z^ad,ad) in the absence 
of data for A^ = 4 (i.e., adiabatic -|- three isocurvature 
modes). In the middle panel the (dashed) curve shows 
the posterior for the prior of section IIB. When this is 
rescaled (by dividing by the posterior with no data and 
then renormalizing) the posterior represented by the solid 
curve results. The right panel shows the results of this 
reweighting on LUf,. With this re- weighting the adiabatic 
fraction increases to Z{ad,ad) = 0.87lg '|'4, with a lower 
isocurvature fraction, /iso = 0.36 ± 0.12. 

V. DISCUSSION 

We have presented a framework within which to con- 
sider correlated mixtures of adiabatic and isocurvature 
perturbations and to obtain simultaneously more model 
independent constraints on cosmological parameters as 
determined from CMB and large-scale structure datasets. 
Applying these methods to recent CMB data, from 
WMAP and several small-scale experiments, indicates 
that current constraints on a single correlated isocurva- 
ture mode are strong (w 10%). These limits, however, de- 
grade substantially when two or three correlated isocur- 
vature modes are allowed, with non-adiabatic fractions 
as large as 60% possible. Including large-scale structure 
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data from the 2dF survey only slightly modifies these 
constraints. 

The parameters significantly modified by including 
correlated isocurvature initial conditions are L0b,T,f3 and 
rig, with the values of the baryon density values twice as 
large as in the pure adiabatic case not ruled out. These 
parameters are strongly correlated with the isocurvature 
amplitudes and for the current data are in fact degenerate 
along a wcU-dcfincd direction in parameter space. Con- 
sequently, the likelihood function is not strongly peaked 
around a given model but relatively flat in a region in- 
terpolating between a set of mixed models and the pure 
adiabatic model. For many of the allowed models the 
large isocurvature fraction did not necessarily accompany 
a substantial decrease in the adiabatic power. This is pos- 
sible because of interference phenomena as we explained. 
While we find that large isocurvature fractions are al- 
lowed we do not find evidence that the inclusion of such 
modes provides a statistically significant better fit to the 
present data. 
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APPENDIX A: RAPID COMPUTATION OF 
COSMOLOGICAL MODELS USING DASH 

1. Motivation and Description 

Cosmological parameter estimation currently requires 
calculating at least 10^ models when a typical 6- 
dimensional parameter space is sampled uniformly. Us- 
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ing an MCMC sampler, which favors the region of high 
hkelihood, reduces the number of computed models to 
Ri 10^ for the same parameter space. Future parame- 
ter estimation studies will aim to measure many more 
parameters as the quality of CMB temperature and po- 
larisation data improves, thereby increasing the compu- 
tational cost. By including non-adiabatic initial condi- 
tions in this study, we have enlarged the parameter space 
to 16 dimensions, requiring computing many more mod- 
els. The development of software that speeds up model 
computation is crucial for sampling such large parameter 
spaces. 

Current codes take between 30 and 60 seconds to com- 
pute a single model, which is prohibitively long for sam- 
pling a large dimensional parameter space. One requires 
a method that both improves the precision of the com- 
puted Ci spectrum as data becomes more accurate with- 
out sacrificing speed. The method adopted here and ex- 
tended is the Davis Anisotropy Shortcut (DASh) [26]. 
With DASh, the accuracy depends primarily on the fine- 
ness of the grid, which depends on the size of available 
fast memory. Consequently, accuracy may be increased 
without sacrificing speed. 

DASh was initially developed as a fast method for pre- 
dicting the angular power spectrum given a set of cosmo- 
logical parameters (the physical baryon density LOb, phys- 
ical matter density uim, relative cosmological constant 
density Qa, relative curvature density CIk and reioniza- 
tion optical depth r) and an arbitrary primordial power 
spectrum. The method relied on the construction of 
grids: a low-^ grid 

for an upper threshold £t['^A, i^x]', a high-£ grid, 

g"[wb,Wm], £>£t; 

a polarization grid 

gi[Wb,Wm]; 

and a reionization grid 

gf[wb,w^,T], 2<e<£R 

for an upper threshold Ir. At each vertex of the low-^ 
and high-^ grids the temperature transfer fmiction Aj{k) 
was pre-computed and stored. Similarly, for the polar- 
ization grid the polarization transfer functions, Af{k), 
were pre-computed and stored. For the reionization grid, 
a function of CJ[t] and CJ[t = 0] was stored at each 
grid point. After the initial construction of the grids, an 
arbitrary C(, spectrum could be calciilatcd in roughly 1 
to 2 seconds, as long as its parameters were within the 
parameter range of the grid. 

The Cf, Cf^, and Cf spectra were calculated by 
interpolating the transfer functions, Aj{k) and Af{k) 
from the values at the nearest vertices on the relevant 



grids. A linear interpolation of the quantities, |A|"(fc)p, 
Aj{k)Af{k) and |A|'(fc)p was performed before inte- 
grating over k to obtain the low Cj spectrum, and the 
high Cj, Cj^ and Cf spectra for r = 0. The high Cj , 
Cj^ and Cf spectra were then rescaled in (. to take into 
account the modified angular diameter distance to recom- 
bination. The low and high £ temperature spectra were 
then combined to obtain the Cj, Cj^ and Cf spectra 
for T — Q. The final spectrum was calculated by mul- 
tiplying by a reionization correction envelope obtained 
from the reionization grid for (. < £r and a suppression 
factor e"^"^ for I > (.r va the temperature case. For the 
temperature-polarization cross-correlation and polariza- 
tion auto-correlation, a fitting function for the reioniza- 
tion bumps was used. More details of the initial release 
[27] of DASh can be found in [26]. Current alternative 
methods to perform rapid CMB model computation or 
to sample parameter space rapidly include [36] and [37]. 

2. Code extensions 

We extended the initial DASh version by including all 
isocurvature modes as in [12]. In the new version, the 
low-^, high-£ and polarization grids were created for each 
pure mode. The cross-correlations between modes were 
calculated at the interpolation stage. We have replaced 
the fitting functions for the reionization bumps with two 
more reionization grids in Cj^ and Cf, for greater ac- 
curacy and because the isocurvature mode bumps had 
very different shapes to the adiabatic mode bumps. The 
three reionization grids were therefore calculated for the 
auto-correlation and cross-correlation modes. 

We have also added the pre-computation of a matter 
transfer function T(fc) grid 

Q'^{k)[uJb,ljJm] 

for each pure mode with = 0.6 for the range in k space 
covered by the 2dFGRS dataset [7]. The matter power 
spectrum was then computed by interpolating T{k) for 
each mode, multiplying by the ratio of growth functions 
for 7^ 0.6 and finally obtaining P{k) by summing over 
modes. The computation of a single P{k) is virtually 
instantaneous (after pre-computation of the grid). 

Finally, due to the 1/k divergences that appear in some 
terms (e.g., the Newtonian potentials) of the line-of-sight 
(LOS) integral [38] for the NIV mode, special care must 
be taken when pre-computing the grids, as this intro- 
duces a numerical error as large as 5% to 10% for < 10 
which could be amplified to as much as 50% after the 
interpolation. We discuss these details next. 

3. Neutrino isocurvature velocity mode 

For the neutrino isocurvature velocity mode, the New- 
tonian potentials which are defined in the conformal New- 
tonian gauge (see for example [39]), diverge as 1/fc and 
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FIG. 21. The different terms in the Line-of-Sight integral (eqn. A5) separated into ones that are positive [left] and negative 
[right] in the limit fc — » 0. The positive terms are : a (solid), &^ (dotted), |§- (dashed) and (long dashed). The negative 
terms are : i) (solid), a (dotted), a (dashed), ^ (long dashed), ^ (dot-dashed), -p- (dot-long dashed) and 11 (short-long 
dashed). Clearly 5^,, f\ and 11 are finite (oc fe) as fe — > whereas the rest diverge as \. 



1/t where t denotes conformal time. To first order they 
are 



* = - 
$ = + 



15 + 4i?,, kt 



1 

TT + ••• 



Ib + ^R^ kt 



1 

TT + . 



(Al) 

(A2) 
(A3) 



where is a constant defined in [12]. These divergences 

have created confusion in the past and have led some 
authors to discard this mode as unphysical. The singu- 
larities above, however, are coordinate singularities, and 
disappear in synchronous gauge. All physical observables 
such as C^s are finite for the NIV mode. 

Let t{) bo the conformal time today, ti an initial con- 
formal time deep into the radiation era, x = k{tQ — t), 
and g{t) = r(t)e'^^*^ the visibility function. The optical 
depth back to t is defined as T{t) = J^* dt f. We use 
the conventions and perturbation variables of [39] ; h and 
7] arc the synchronous metric perturbations, S-y the pho- 
ton density contrast, Oh the baryon velocity divergence, 
n = -l- Go -h G2 the polarization source term with 
and Ge the sum and difference of the two photon polar- 
izations, and a = -^{h + 6?)) the gauge transformation 
variable. For the neutrino isocurvature velocity mode, at 
early times and on large scales the above variables are 
given to leading order in t by 

5-^ (xkt + . . . , 9b oc k + . . . , T] oc kt + ... , 



1 

a oc — + 

k 



U(xkt+.... (A4) 



The line-of-sight integral [38] for the temperature 
transfer function is given by 



Mk) = J^ dtje{x)ie^{ri + a) 



+ 9{t) 



Ob 3 



16fc2 



n 



(A5) 



Most terms in the above equation diverge as 1/k as 
shown in Fig. 21. The transfer functions, Af:{k), how- 
ever, are finite, since the Bessel functions are well ap- 
proximated by je{kt) oc {ktY at small k, which is suffi- 
cient to cancel the 1/fc divergence. This means that if 
A^(A:) is calculated using an exact Boltzmann code, it 
will be finite. Using the LOS integral above, however, 
one should worry that in computing Ae{k) the divergent 
terms will cause numerical noise on large scales, because 
of the subtraction of large numbers of (nearly) equal mag- 
nitude. This indeed occurs, with resulting errors ranging 
from 5% to 10%. Such errors are less a problem for codes 
based entirely on the LOS method because the errors lie 
below the cosmic variance. Using the DASh code, how- 
ever, the interpolation between models with this type of 
random error boosts the final error to as much as 30% to 
50% at £ < 10 in some cases. 

In order to eliminate the above numerical error, one 
can integrate eqn. (A5) by parts, so that 



Ae{k) 



dt yeix) 

dji 
dx 



1 



e^rj + -g{t) [5^ + -^ 



1, 



e^M+^(0^n 
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FIG. 22. Accuracy of the computation of the temperature transfer function. [Loft] The temperature transfer 
function A2(fe) calculated using a Boltzmann code (solid), with the line-of-sight integral (dashed) and with the integration 

by parts (dotted). [Upper Right] The angular power spectrum Ce for a Boltzmann code (solid), LOS (dashed) and by parts 
(dotted). [Lower Right] Percent difference between LOS and Boltzmann (solid) and by parts and Boltzmann (dashed). 



4. Accuracy and Timing 




FIG. 23. The temperature transfer functions A2oo(fc) cal- 
culated using the LOS method (dotted) and linear interpola- 
tion (solid). The prominent feature is a beat of much lower 
frequency than the transfer function. This is what limits 
the interpolation accuracy. (Note that the actual code in 
DASli interpolates between the squares of the transfer func- 
tions which decreases the effect of the beat). 



+ g(t) [ka+-f + 



-n 



k 16fc 



A quick check demonstrates that all the terms in the 
above equation are finite as k 0, which gives stable 
and accurate results when integrated numerically. We 
compare the different integration methods in Fig. 22. 



The LOS method is an extremely accurate method and 
is limited only by the accuracy of the Besscl functions, 
which for flat models is 0.1%, as was also demonstratcid in 
[40]. Indeed this is the accuracy of all the transfer func- 
tions computed in DASh. The approximation techniques 
introduce further errors, namely interpolation error for 
all models and approximation error on large scales for 
reionized models due to the approximate method used to 
include reionization effects on large scales. 

The interpolation error is due to an intrinsic 0.1% error 
of the transfer functions which is amplified by interpola- 
tion, and a beat phenomenon which arises because the 
frequencies of the transfer functions being interpolated 
are slightly offset from each other and the target function. 
The intrinsic interpolation error could be remedied if the 
transfer functions were calculated with an exact Boltz- 
mann code, though this would be at a cost in speed of pre- 
computation. To illustrate the effect of the interpolation 
error we plot the calculated LOS transfer function and 
the interpolated one in Fig. 23. Note that the interpola- 
tion in Fig. 23 is between transfer functions and not their 
squares, as is actually done in DASh, to dramatize the 
effect. Interpolating the squares of the transfer functions 
reduces this error, at the expense of having to interpolate 
3 X 10 transfer functions (counting the temperature and 
polarization for all auto- and cross-correlation modes) 
instead of 2 x 4 (where the cross-correlation modes are 
calculated after interpolation). 

To check that the inclusion of additional initial con- 
ditions, and the extension to compute the matter power 
spectrum, were sufficiently accurate wc tested the DASh 
code against 2048 models sampled from a grid in 
{wbjWcf^AjT}, in the case of CMB temperature and 
polarization, and 512 models sampled from a grid in 
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FIG. 24. Maximum and RMS percentage differences between the DASh interpolated spectra and the 

line-of-sight computed spectra. Maximum (dashed) and RMS (sohd) percentage differences of DASh compared to the 
Line-Of-Sight method. From left to right we show the CMB temperature spectrum, CMB polaxization spectrum and matter 
power spectrum P{k). From top to bottom we show the adiabatic, CDM isocurvature, neutrino isocurvature velocity and 
neutrino isocurvature density modes. The errors for Ce were calculated using 2048 models and for P{k) using 512 models (we 
excluded variation of t for the matter power spectrum. 
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{ujh,LUc,^A}, in the case of the matter power spectrum. 
We plot the maximum and root-mcan-square (rms) er- 
rors in Fig. 24 for all auto-correlation mode spectra. For 
the CMB temperature spectrum (for which the most ac- 
curate measurements exist) we see that the maximum 
error is loss than 2% over a wide range in £, which is bet- 
ter than the WMAP measurement error. On large scales 
{£ < 10) the maximum error increases to 10% for some 
modes but this is still well below cosmic variance. 

We tested the effect of these numerical errors on the 
likelihood computed for each model and found the mea- 
sured difference to be negligibe. 

The isocurvature-modified DASh described here took 
approximately 3-4 seconds on a IGHz Pentiumlll pro- 
cessor to compute the CMB temperature, polarization, 
and cross-correlation spectra as well as the matter power 
spectrum, for a single mode. The full matrix of 4 modes 
and their cross-correlations required approximately 15 
seconds per model. 
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